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The zero-temperature linear response conductance through an interacting mesoscopic region at- 
tached to noninteracting leads is investigated. We present a set of formulas expressing the con- 
ductance in terms of the ground-state energy of an auxiliary system, namely a ring threaded by a 
magnetic flux and containing the correlated electron region. We prove that the formalism is exact 
if the ground state of the system is a Fermi liquid. We show that in such systems the ground-state 
energy is a universal function of the magnetic flux, where the conductance is the relevant parameter. 
The method is illustrated with results for the transport through an interacting quantum dot and a 
simple Aharonov-Bohm ring with Kondo-Fano resonance physics. 



I. INTRODUCTION 

In the last decade technological advances enabled con- 
trolled fabrication of small regions connected to leads and 
the conductance, relating the current through such a sys- 
tem to the voltage applied between the leads, proved to 
be the most important property of such systems. There 
is a number of such examples, e.g. metallic islands pre- 
pared by e-beam lithography or small metallic grains, * 
semiconductor quantum dots^ or a single large molecule 
such as a carbon nanotube or DNA. It is possible to 
break a metallic contact and measure the transport prop- 
erties of an atomic-size bridge that forms in the break^ 
or even measure the conductance of a single hydrogen 
molecule. 4 Recent measurements of conductance through 
single molecules proved that strong electron correlations 
can play an important role in such systems £ 

The transport in noninteracting mesoscopic systems 
is theoretically well described in the framework of the 
Landauer-Biittiker formalism. The conductance G is at 
zero temperature determined with the Landauer-Biittiker 
formula^ 



G = G \t(e F )\ 2 , G o-— • (1) 

The key quantity here is the single particle transmission 
amplitude t(ep) for electrons at the Fermi energy. The 
formula proved to be very useful and reliable, as long 
as electron-electron interaction in a sample is negligible. 
However, the Landauer-Biittiker formalism cannot be di- 
rectly applied to systems where the interaction between 
electrons plays an important role. Several approaches 
have been developed to allow one to treat such systems. 
The Kubo formalism provides us with a conductance for- 
mula which is applicable in the linear response regime 
and was intensively studied by Oguriii& A more general 
approach applicable also to non-equilibrium cases was 
developed by Meir and Wingreen^ Recently, ab initio 
methods to study the transport through small molecular 
junctions were also applied^ 



II. CONDUCTANCE FORMULAS FOR FERMI 
LIQUID SYSTEMS 

The relevant system is schematically presented in 
Fig. nja). A mesoscopic interacting region, which could 
be a molecule, a quantum dot, a quantum dot array or a 
similar 'artificial molecule' system J^s attached to nonin- 
teracting leads. As shown in Ref. O (hereafter referred 
to as RR) , the conductance of such a system can be de- 
termined solely from the ground-state energy of an aux- 
iliary system, formed by connecting the leads of the orig- 
inal system into a ring and threaded by a magnetic flux, 
Fig. IHb). The main advantage of this method is the 
fact that it is often much easier to calculate the ground- 
state energy (for example, using variational or quan- 
tum Monte Carlo methods) than the Green's function, 
which is needed in the Kubo and Keldysh approaches. 
The method is applicable only to a certain class of sys- 
tems, namely to those exhibiting Fermi liquid properties, 
at zero temperature and in the linear response regime. 
However, in this quite restrictive domain of validity, the 
method promises to be easier to use than the methods 
mentioned above. 

The basic property that characterizes Fermi liquid sys- 
tems is that the states of a noninteracting system of 
electrons are continuously transformed into states of the 
interacting system as the interaction strength increases 
from zero to its actual valued One can then study the 
properties of such a system by means of the perturba- 
tion theory, regarding the interaction strength as the per- 
turbation parameter. Dynamics of Fermi liquid systems 
at low temperature and in the linear response regime 
is governed by quasiparticles. However, the question 
how quasiparticles propagate in a correlated system is 
a non-trivial one. The answer can be extracted from the 
Green's function for a particular problem if it is known. 
An alternative way, which we advocate in this paper, is 
to analyze the excitation spectrum of a system directly. 
If Em and Em+i are the ground-state energies of an in- 
teracting ring system containing M and M + 1 electrons, 
respectively, the energy difference can be attributed to 
the first quasiparticle energy level e above the Fermi en- 
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Figure 1: (a) Schematic picture of a sample with interaction 
connected to noninteracting leads, (b) The sample embedded 
in a ring formed by joining the left and right leads of the 
system (a). Auxiliary magnetic flux 3> = -(f) penetrates the 
ring. 

ergy, 

e = Em+i — Em- (2) 

The variation of the quasiparticle energy with flux 
threading the ring allows us to determine the conduc- 
tance of the system. The complete proof of the formal- 
ism is given in RR and a brief overview is presented in 
the next Section. Here we show how the method can be 
implemented in practice. 

The key property of ring systems presented in Fig.njb) 
is the universality expressed in the variation of the 
ground-state energy with auxiliary magnetic flux through 
the ring. Here we assume a system obeys the time rever- 
sal symmetry. The more general case is presented in the 
last Section. For an even number of electrons in the sys- 
tem and a large number of sites in the ring N — > oo the 
ground-state energy takes a universal form 

E(4>)-E(^j = ^(arcco S 2 (TV9cos0)-^), (3) 

where the average level spacing at the Fermi energy A = 
[Np (ei?)] _1 is determined by the density of states at the 
Fermi energy in an infinite noninteracting lead p (ep) and 
g = G/Go is the dimensionless conductance. For systems 
with an odd number of electrons, the ground-state energy 
is given with 

(|) = ^arcsin 2 (Vgcos0). (4) 

It should be mentioned that the ground-state energy of 
an interacting ring system exactly corresponds to the ex- 
pression for persistent currents in noninteracting rings, 



as derived by Gogolin^ The only parameter determin- 
ing the ground-state energy is the conductance g of the 
original system, Fig. Ufa). I n Fig. the ground-state en- 
ergy as a function of the flux <j> for even and odd numbers 
of electrons is presented. It should be pointed out that 
as successive quasiparticle levels are being occupied, the 
points = and (/) = ir interchange their roles (as is also 
the case for noninteracting systems), and that the peri- 
odicity in even and odd cases are ir and 7r/2, respectively. 
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Figure 2: The ground-state energy of an interacting system as 
a function of flux (j> for an even (a) and an odd (b) number of 
electrons for g going from to 1 in steps of 0.1 and N — > oo. 

If the ground state of the system in Fig.QIb) is known, 
the conductance of the original open system can be ex- 
tracted from Eq. (3J) [or Eq. J3J]. There are several ways 
how to determine g from Eq. <EJ. The simplest seems to 
be the use of the relation 

, fir EM -E(0)\ 
g = Sm U A )' (5) 

where E(0) and E(ir) are the ground-state energies of the 
ring system for cj) — and 4 1 — n, respectively. The first 
advantage of this formula is the fact that the energies 
can be calculated using periodic and antiperiodic bound- 
ary conditions, respectively, and thus the wave functions 
of the system can be taken real. Additional advantage 
is fast convergence with N, as briefly discussed in the 
last Section. This formula was derived in the g — > 
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limit as g = ( ^ [E (tt) - E (0)]) by Favand and Mila for 
noninteracting systems and applied to interacting Hub- 
bard chains^ More recently, a similar approach was per- 
formed in Ref . ITH In Fig.OJa) the use of formula Eq. © 
is schematically presented. 
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Figure 4: Conductance vs. charge stiffness using D at the 
energy minimum (full line) and D at the maximum (dashed 
line). Note the quadratic dependence g oc D 2 for g — > 0. 



in the minimum of the energy vs. flux curve, D = 
Y d 2 E/d(j) 2 \ E _ . |21 The sensitivity of the ground state 
energy to flux has been applied also in the context of 
electron localisation. 22 One can also define the corre- 
sponding quantity for the energy maximum as D = 
— d 2 E/d<f> 2 \ . From Eq. (0 the conductance can 

be related to the charge stiffness with an implicit relation, 



1 d 2 E 
A ~W 



±- 



E— min, max 



1-5 



arccos (±y/g) . (7) 



Here the upper and the lower signs correspond to the 
second derivative at a minimum and at a maximum of 
the energy vs. flux curve, respectively. In general, this 
equation has to be solved numerically and the solutions 
are presented in Fig. Q] In the limit of a very small con- 
ductance and in the vicinity of the unitary limit, analytic 
formulas are available 



Figure 3: The conductance can be extracted from: (a) the 
two-point formula, Eq. 10, (b) the persistent current formula, 
Eq. 10, or (c) the charge stiffness formula, Eq. 0. 

The derivative of the ground-state energy with respect 
to flux gives the persistent current in the ring j(</>) = 
j-^^L The second formula relates the conductance to 
the persistent current at 



2 > 



7T H 

9 = 1 Ae J I 2 



(6) 



This relation was recently derived for noninteracting 
systems^ and successfully applied to systems with 
interactioniA2iA£ 

At T = the charge stiffness is an important quantity 
describing the charge transport in correlated systems 31 
It is defined as the second derivative of the ground- 
state energy of the system with respect to the flux 
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[2vp{e F )D] 2 , g 
(l + ^[27rp(e F )D]) 2 , g 



(8) 



Note that there is a quadratic relation between the con- 
ductance and the charge stiffness in the low conductance 
limit. The corresponding formulas for the maximum of 
the energy vs. flux curve are 



2np {e F )D 



.9^0, 



1 - 



fJ 



1. 



(9) 



It should be stressed that the validity of all formulas 
presented in this Section is based on an assumption that 
the number of sites in the ring is sufficiently large accord- 
ing to the condition^ 



1 dy/g(e F ) 
p{e F ) de F 



(10) 
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This means that if g (ep) exhibits sharp resonances, as is 
the case, e.g., in chaotic systems^ the calculation has to 
be performed on a large auxiliary ring system and in such 
cases the method might be impractical compared to other 
methods. On the other hand, for systems with strong 
interaction the method promises to be extremely efficient 
already for ring systems of a moderate size i 1 ^ 1 ^ 24 



III. PROOF OF THE FORMALISM 

The complete proof of the formalism is presented in 
RR, here we briefly describe the main steps. The proof 
strongly relies on an assumption that the ground state of 
the system under investigation is a Fermi liquid^ 

We start with the linear response conductance of a gen- 
eral interacting system of the type shown in Fig. da). 
The conductance can be calculated from the Kubo for- 
mula^ 



ITT 



lim 

w^o to + id 



:fl 



i5). 



(11) 



where Uji (lo + iS) is the retarded current-current corre- 
lation function. For Fermi liquid systems, the current- 
current correlation function can be calculated within the 
perturbation theory. At T — 0, only the bubble diagram 
gives a non-vanishing contribution^ and the conductance 
can be expressed in terms of the Green's function G n > n (z) 
of the system, 



9 = 



-mp(e F ) 



-ik F (n'-n) Gntn {ep+iS) 



(12) 



where n and n' are sites in the left and the right lead, 
respectively. 

In Fermi liquid systems obeying the time-reversal 
symmetry^ the imaginary part of the retarded self- 
energy at T — vanishes at the Fermi energy and is 
quadratic for frequencies close to the Fermi energyiSl Us- 
ing the Fermi energy as the origin of the energy scale, i.e. 
lo — ep — * to, we can express this as 



ImSy (lo + iS) oc lo 2 



(13) 



Close to the Fermi energy, the self-energy can be ex- 
panded in powers of lo resulting in an approximation to 
the Green's function, 



G 1 (lo + id) = loI - H (0) - S (0 + id) 



as (lo + is) 



duo 



O(oo 2 ) 



(14) 



Here H^ ^ contains matrix elements of the noninteracting 
part of the Hamiltonian. The Green's function for lo close 
to the Fermi energy can then be expressed as 



G 1 (lo + iS) = Z-^GT 1 (lo + iS) Z" 1 / 2 + O (lo 2 



(15) 



where we defined the quasiparticle Green's function 

G 1 (lo + iS) = loI - H (16) 

as the Green's function of a noninteracting quasiparticle 
Hamiltonian 



H = Z 1 / 2 



H(°! + S (0 + id) 



z 1 / 2 , 



(17) 



and introduced the renormalization factor matrix Z. Ma- 
trix elements of Z differ from those of an identity matrix 
only if they correspond to sites of the central region. 

The reason for introducing the quasiparticle Hamilto- 
nian is to obtain an alternative expression for the con- 
ductance in terms of the quasiparticle Green's function. 
Eq. Ill ."it relates the values of the true and the quasipar- 
ticle Green's function at the Fermi energy, 



G (0 + iS) = Z 1/2 G (0 + iS) Z 1/2 . 



(18) 



Specifically, G n > n (0 + iS) = G n i n (0 + id) if both n and 
n' are sites in the leads, as a consequence of the prop- 
erties of the renormalization factor matrix Z discussed 
above. Eq. lfl2ll then tells us that the zero-temperature 
conductance of a Fermi liquid system is identical to the 
zero-temperature conductance of a noninteracting system 
defined with the quasiparticle Hamiltonian for a given 
value of the Fermi energy. 

These conclusions are valid if the central region is cou- 
pled to semi-infinite leads. Here we generalize the con- 
cept of quasiparticles to a finite ring system with N sites 
and M electrons, threaded by a magnetic flux 4>. One can 
define the quasiparticle Hamiltonian for such a system, 



H (N, 0; M) = Z 1/2 [H (0) (N, 4>) + S (0 + iS) 



7}' 2 . 



(19) 

Here the self-energy and the renormalization factor ma- 
trix are determined in the thermodynamic limit where, 
as we prove in RR, they are independent of <j> and corre- 
spond to those of an infinite two-lead system. 

Suppose now that we knew the exact values of the 
renormalized matrix elements in the quasiparticle Hamil- 
tonian Q19J1 , As this is a noninteracting Hamiltonian, we 
could then apply the conductance formulas presented in 
the previous Section (the proof of validity of energy for- 
mulas for noninteracting systems is given in RR and the 
corresponding result for persistent currents in Ref. . 
Such a procedure would provide us with the exact con- 
ductance of the original interacting system. However, to 
obtain the values of the renormalized matrix elements, 
one needs to calculate the self-energy of the system, 
which is a difficult many-body problem. In RR we study 
the excitation spectrum of a finite ring system with in- 
teraction and threaded with a magnetic flux. We show 
that 



E (N, 4>;M + 1) — E (N, 0; M) = 
= e(iY,^M;l)+e>(iV-i) , 



(20) 
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where E (N, cj>\ M) and E (N, 0; M + 1) are the ground- 
state energies of the interacting Hamiltonian for a ring 
system with ./V sites and flux (j>, containing M and M + 1 
electrons, respectively, and e (N, </>; M; 1) is the energy of 
the first single-electron level above the Fermi energy of 
the finite ring quasiparticle Hamiltonian iflfljl . The error 
in Eq. l|20j) is small enough that the proof of the ground- 
state energy formulas for noninteracting systems, which 
involves only the properties of a set of neighboring single- 
electron energy levels, remains valid also for interacting 
Fermi liquid systems, provided a system is a Fermi liquid 
for all values of the Fermi energy below its actual value. 



IV. EXAMPLES 



A. Noninteracting system 
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In this Section we discuss the convergence properties 
of the conductance formulas. As the first example we 
take a double-barrier potential scattering problem pre- 
sented in Fig.0 Results of various formulas for different 
number of sites in the ring are presented in Fig. U3 The 
exact zero-temperature conductance for this system ex- 
hibits a sharp resonance peak superimposed on a smooth 
background conductance. We notice immediately that as 
the number of sites in the ring increases, the convergence 
is generally faster in the region where the conductance 
is smooth than in the resonance region, which is consis- 
tent with the condition Eq. Comparing the results 
obtained employing the two-point formula Eq. Q and 
the persistent current formula Eq. JSJ we observe that 
the convergence is better for the two-point formula ex- 
pressing the conductance in terms of the difference of the 
energies at = and cj) — n. From the computational 
point of view there is an additional advantage of the two- 
point formula. In this case, all the matrix elements can 
be made real if one chooses such a vector potential that 
only one hopping matrix element if modified by the flux 
as then the additional phase factor is e ±l7r = — 1. 




Figure 5: A double barrier noninteracting system. The height 
of the barriers is 0.5t, where t is the hopping matrix element 
between neighboring sites. 



Figure 6: (a) The conductance through the system in Fig. [£] 
calculated from the two-point formula, Eq. 10, and (b) from 
the persistent current formula, Eq. J§J|. Note the different 
convergence behavior of the two formulas. 

are compared to exact conductance of the Bethe ansatz 
approach^ To calculate the conductance, Eq. (0 was 
used, with the ground-state energies at <j> — and <j> = tt 
obtained using a variational method described in RR. For 
each position of the Cd level relative to the Fermi energy, 
we increased the number of sites in the ring until the 
conductance converged. The number of sites needed to 
achieve the convergence was the lowest in the empty or- 
bital regime and the highest (about 1000 for the system 
shown in Fig. |HJ in the Kondo regime. This is a conse- 
quence of Eq. ifTfljl as a narrow resonance related to the 
Kondo resonance appears in the transmission probabil- 
ity of the quasiparticle Hamiltonian Ill7fl in the Kondo 
regime. 




Figure 7: The Anderson impurity model realized as a quan- 
tum dot coupled to two leads. The dot is described with the 
energy level td and the Coulomb energy of a doubly occupied 
level U . t' is the hopping between the dot and leads. 



C. Interacting Aharonov Bohm rings 



B. Anderson impurity model 

As a nontrivial example of the use of the formalism 
we calculate the zero-temperature conductance of a sin- 
gle impurity Anderson model realized as a quantum dot 
attached to leads as shown in Fig. [3 In Fig.|H|the results 



One can generalize the conductance formulas to sys- 
tems which exhibit time reversal asymmetry, such as 
is e.g. an Aharonov-Bohm (AB) type of system pre- 
sented in Fig. |2fa)i^ If there is no AB flux thread- 
ing the mesoscopic region, the time reversal symme- 
try is restored and the energy is an even function 
of <p. In the general case, the energy extremum is 
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Figure 8: The zero-temperature conductance calculated from 
ground-state energy vs. magnetic flux in a finite ring system 
using the variational method described in RR with 3 and 7 
basis functions. For comparison, the exact Bethe ansatz result 
is presented with a dashed line. The system is shown in Fig.Q 
with U = 0.64t and t' = 0.2t. 



shifted to a non-trivial point 0o (^f), as illustrated in 
Fig. injc) . The ground-state energy is then generalized 
to E(4>) = 7r~ 2 A arccos 2 (Ty^ cos ~ ( £ f)]) + const, 
for an even number of electrons in a system and to 
E((f>) — 7r _2 A arcsin 2 (y^cos [0 — 0o (ef)]) + const, for 
an odd number of electrons in a system. From the for- 
mer expression, the transmission probability can be ex- 
tracted, and the conductance is given by^ 



.9 



sin — - 
2 



ttE{ 



n)-Ei 



(21) 



where (f>o = (f>o (ep) is determined by the position of the 
minimum (or maximum) in the energy vs. flux curve. 
The conductance can also be calculated from the more 
convenient four-point formula 32 



.9 = sm 



2 fn E(tt)-E(0) 
A 



-sin 



7T E{n/2)-E{-ir/2) 



A 



and 0o is determined with the expression 



= — arctan 



3in(f 



E(ir/2)-E(—n/2) 
A 



sin( f MZL)_M0)) 



(22) 



(23) 



If no AB flux is present in the mesoscopic region, we 
recover the two-point formula Eq. © since E(-k/2) = 
E (—it/2) in this case. 

If the time reversal symmetry is broken due to AB flux, 
Eq. ijlSj) is not valid and the proof of the previous Section 
has to be reconsidered. Repeating the steps as presented 
in detail in RR, the proof is restored and basically un- 
changed if the self energy obeys the relation 



1 

2? 



[Ejj (uj + iS) — T,ij (uo — i8)\ oc ui 2 



(24) 



It follows that the linear response conductance of an in- 
teracting AB system at zero temperature is given by the 




Figure 9: (a) Interacting mesoscopic region (gray-shaded 
sites), threaded by magnetic flux and coupled to noninteract- 
ing leads, (b) Auxiliary ring system, (c) Behavior of energy 
levels as the flux threading the ring is varied, (d) Example 
system with a quantum dot embedded in an AB ring. 



four-point formula Eq. (22j|. The condition Eq. l|24ll is 
fulfilled if the system is a Fermi liquid. 

In order to demonstrate the practical value of the 
method, we quantitatively analyze the conductance 
through an Aharonov-Bohm ring with a quantum dot 
placed in one of the arms^fliSI as presented in Fig. Efcl) . 
In Fig.E3 a convergence test of the method is shown. The 
convergence with N is fast in the empty orbital regime 
and becomes progressively slower as shifts toward the 
Kondo regime, as was also the case for a quantum dot at- 
tached to leads. The converged conductance curve is in 
excellent agreement with the numerical renormalization 
group result of Ref . El 



V. SUMMARY 

We have demonstrated how the zero-temperature con- 
ductance of a sample with electron-electron interaction, 
attached to noninteracting leads can be determined. The 
method is extremely simple and is based on several for- 
mulas relating the conductance to the ground-state en- 
ergy of an auxiliary ring system. The conductance is 
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Figure 10: Zero-temperature linear response conductance of 
the system in Fig. Eld) as a function of level position Ed for 
various number of ring sites N and 4>ab = 7r/4. The dotted 
line is the NRG result from Ref.lsil Parameters: t\ = 0.177to 
(r = 0.125t ), ti = 0.298to, U = t = 8I\ 

determined from the ground-state energy of an interact- 
ing system, while in more traditional approaches, one 
needs to know the Green's function of the system. The 
advantage of the present method comes from the fact 
that the ground-state energy is often relatively simple to 
obtain by various numerical approaches, including varia- 



tional methods. Let us summarize the key points of the 
method: 

(1) The "open" problem of the conductance through a 
sample coupled to semi-infinite leads is mapped on to a 
"closed" problem, namely a ring threaded by a magnetic 
flux and containing the same correlated electron region. 

(2) For the case of a Fermi liquid interacting system, 
even with broken time-reversal symmetry, it is shown 
that the zero-temperature conductance can be deduced 
from the variation of the ground state energy with the 
flux in a large, but finite ring system. 

(3) In order to prove this, the concept of Fermi liquid 
quasiparticles is extended to finite, but large systems. 
The conductance formulas give the conductance of a sys- 
tem of noninteracting quasiparticles, which is equal to 
the conductance of the original interacting system. 

(4) The results of our method are compared to results 
of other approaches for problems such as the transport 
through a quantum dot containing interacting electrons. 
The comparison shows an excellent quantitative agree- 
ment with exact Bethe ansatz results. We have demon- 
strated the usefulness of the formula also by applying it 
to a prototype system exhibiting Kondo-Fano behavior. 
Results based on the four-point formula confirm the re- 
sults of the numerical renormalization group method. 
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